
HEAT AND MASS TRANSFER IN BINARY INERT GAS FLOW FOR 
DISTRIBUTIONS OF TEMPERATURE AND CONCENTRATION 
RENDERING THE PROPERTIES NEARLY CONSTANT 


Ernst W. Adams, John D. Warmbrod, C. Lee Fox, and Robert M. Huffaker 


/ '/ • 



MSFC - Form 523 (Rev. November I960) 


GEORGE C. MARSHALL SPACE FLIGHT CENTER 


MTP-AERO-64- 2 


HEAT AND MASS TRANSFER IN BINARY INERT GAS FLOW FOR 
DISTRIBUTIONS OF TEMPERATURE AND CONCENTRATION 
RENDERING THE PROPERTIES NEARLY CONSTANT 

By 

Ernst W. Adams, John D. Warmbrod, C. Lee Fox, 
and Robert M. Huffaker 


ABSTRACT ^ yv 

/-i 

Injection of a foreign gas into the laminar air boundary layer is 
considered. The mixture properties are arbitrary functions of tempera- 
ture T and foreign gas concentration w. Unless the properties are con- 
stant, similarity transformations are valid only at the stagnation point 
or for a wall at constant pressure. Solutions of the three similarity 
equations are quite cumbersome because triple iterations are involved 
at the wall to satisfy three conditions at the outer edge of the boundary 
layer. A correlation formula of rigorous numerical solutions for q w /q wo 
has been derived in Reference 14 for the constant pressure case. Here, 
q w is the heat transfer in the presence of mass transfer, whose absence 
is denoted by subscript o. This correlation formula is linear in the 
mass transfer rate p w v w , and, therefore, fails for large values of p w v w . 

A simple engineering solution method of the similarity equations 
is proposed here, which gives those points of the q w /q wo versus p^v^^ 
relationship for which five property parameters are nearly constant 
across the boundary layer. If these five conditions are satisfied, 
the differential equations in similarity variables can be uncoupled. 

By use of an auxiliary graph, the momentum equation can be integrated 
directly as an initial value problem. This solution is used to deter- 
mine q w by quadratures. Correlation formulas for q w are presented for 
both the constant pressure and the stagnation point cases. 

The five conditions on the mixture properties can be satisfied in 
an approximate way for the injection of H 2 0, H e , or H 2 into air, provided 
dissipation effects are sufficiently small. Comparison to the correla- 
tion formula of Reference 14 shows very good agreement if the Mach number 
M^ = 0, and some difference for Moo = 3. The results for H 2 0-air mixtures 
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cover the range of q w /q wo values from unity to values as low as 0.5. 
Only very small injection rates are compatible with the five conditions 
if H e or H 2 is injected into air. The theory is worked out in this 
paper for the case of constant pressure and air as the primary flowing 
mediu “- 
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I. INTRODUCTION 


The problem of simultaneous heat and mass transfer began to attract 
considerable attention about ten years ago. A two-way interaction between 
these transfer rates occurs in case of ablating surfaces, and a one- 
directional interaction in case of transpiration cooling. Ablation is 
experienced by objects traversing an atmosphere at sufficiently high 
speed; ablation also takes place under suitable conditions if a medium 
flowing through a chemical reactor undergoes a reaction with a solid or 
a liquid deposited in the reactor. Transpiration cooling of the surfaces 
of wings, gas turbine blades, etc., has been studied extensively in the 
past. 


Even though the assumptions upon which boundary layer theory is 
based are violated in case of normal mass transfer across a wall at rest, 
experiments confirm pertinent boundary layer solutions within measuring 
accuracy. Laminar boundary layer solutions are physically valid only 
for sufficiently small injection rates normal to the surface because of 
stability considerations. Solutions of the stationary boundary layer 
equations may be used to approximate nonstationary flight situations. 

An exceptional case occurs when the ablation rate depends not only on 
the instantaneous heat transfer rate but also on the transient conduc- 
tion of heat in the interior of the wall whose surface is ablating (see 
Reference 1). If time-dependent effects may be ignored, the double inter- 
action between heat and mass transfer can be analyzed for a given flight 
state by interpolating between a number of simple- interaction solutions 
with prescribed mass injection. 
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Theoretical contributions to the heat and mass transfer problem under 
discussion deal predominantly with stationary, binary, and laminar boundary 
layer flow consisting of a main ("primary") stream and an injected ("foreign") 
constituent. This problem is governed by a system of four simultaneous dif- 
ferential equations, (1) through (4) in Section 2, representing conservation 
of matter, momentum, foreign species, and energy. Up to now, these equa- 
tions have been solved rigorously only under simplifying assumptions on the 
mixture properties and/or the flow parameters. Inspection of equations (1) 
through (4) yields the conditions under which either the mass concentration 
w(x, y) of the foreign gas is determined as a function of the temperature 
distribution T(x, y) or the distributions of both w(x, y) (see Reference 3, 
p. 608 and Reference 9, p. 458). Similarity solutions of equations (1) 
through (4) have been published for both constant and variable mixture pro- 
perties. If the Mach number, the temperature difference across the boundary 
layer, and the foreign gas concentration are sufficiently small, constant 
mixture properties may be assumed. This case has been treated in References 
15 and 25. Similarity solutions then hold true for any Euler number m. If 
the mixture properties are allowed to depend on temperature and foreign gas 
concentration, similarity solutions exist only for the stagnation point 
(m = 1) and for zero pressure gradient (m = 0). The case m = 1 is dealt 
with in References 16 and 17; the case M = 0 is treated in References 11 
and 14. In some references, homogeneous or heterogeneous chemical reactions 
of assumed order between the constituents are included in the discussion, 
e.g,. Chapter 19 of Reference 3. Solutions for more than two chemical con- 
stituents have been derived in connection with dissociation or ionization of 
the constituents in the mixture, e.g., References 21 and 23. 

Even if a similarity transformation has been introduced, the system 
of simultaneous nonlinear differential equations offers formidable resis- 
tance to numerical solutions, provided the mixture properties are allowed 
to depend on both foreign gas concentration w and temperature T as is pre- 
dicted by the kinetic theory of gases, e.g.. References 6 and 19. This 
situation becomes evident by considering that three boundary conditions 
at the outer edge of the boundary layer have to be satisfied by simultaneous 
iterations for the gradients of velocity, temperature, and concentration at 
the wall. The momentum equation in similarity variables becomes independent 
of the differential equations for T and w if the produce of viscosity p and 
density p is assumed to be constant. However, for the compressible boundary 
layer flow of a single gas, similarity solutions show a significant difference 
if the case of constant pp is dependent on T according to the equation of 
state and Sutherland 1 s viscosity law for air, e.g.. Chapter 4 of Reference 27. 

It is shown in this report for mixtures consisting of air and H 2 0 
vapor, He, or H 2 that it is possible for any given set of boundary values 
T w and Too to determine the wall concentration w w so that the relevant 
mixture properties take approximately the same values at the wall and at 
the outer edge of the boundary layer. Under suitable conditions specified 
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in this report, these mixture properties then must be nearly constant 
across the entire boundary layer. If it may be assumed that pp, = const, 
in the boundary layer, the momentum equation in, similarity variables can 
be solved independently of the other differential equations. The equa- 
tion for the foreign gas concentration and, subsequently, the one for 
the temperature distribution, can then be solved by quadratures involving 
the solution of the momentum equation. Evidently, the numerical work of 
this solution method is very much smaller than the work needed to inte- 
grate the simultaneous system of three similarity differential equations 
with split boundary conditions. This is particularly true since experi- 
ence shows that numerical instabilities in these integration procedures 
limit their range of applicability. Also, it is impossible to check the 
accuracy of numerical integrations of the differential equations by means 
as simple as the ones needed to check the quadratures involved in the 
method proposed in this paper. 


II. THE DIFFERENTIAL EQUATIONS AND BOUNDARY CONDITIONS 

According to page 457 of Reference 9, the stationary, two-dimensional, 
and laminar boundary layer flow of a binary, compressible mixture is 
described by the following system of differential equations if both pres- 
sure gradient and thermal diffusion are negligible and if the mixture is 
chemically inert: 


d(P u ) + d ( pv) = o 

Bx + By 


pu s + w If * h G 37) ■ °- (2) 

pU lx + pv 3y " Jy ( D |y) " °’ and <3) 



C_v ^2: B 

p p By " By 




PD < c p 
*a 



Bt Bw 

3737 ’ 0 


( 4 ) 


In the equations (2), (3), and (4), the third term represents, respectively, 
the diffusion of momentum, matter, and heat. The last two terms in equa- 
tion (4) stand for dissipated mechanical energy and for internal energy 
exchanged due to the diffusion flows of primary and injected medium, 
respectively. 
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The following boundary conditions are specified: 

u(x, 0) = 0, v(x, 0) = v w (x), w(x, 0) = w w (x), and T(x, 0) = T w (x) 

(5) 

at the wall y = 0; and 


lim u = Uoo(x), 11m w = Woo(x) 5 and lim T = Too(x) ... 

y — > oo y — ^ co y — > oo 

at the outer edge of the boundary layer. The quantities v w , w w , u^, w^, 
and Too have to be given. Since the wall is supposedly impermeable to 
the primary gas, convection of this gas toward the wall is balanced by 
diffusion of this gas in the opposite direction, i.e., 


Pw v w (! ~ w w ) - ~ Pw°w ( J • ( 7 ) 

v 'w 

Provided suitable initial conditions have been prescribed at some upstream 
station x = const., e.g., at the stagnation point x = 0, equations (1) 
through (6) completely determine the functions u(x, y) , v(x, y) , w(x, y) , 
and T(x, y) . In the following, this system is solved for the tempera- 
ture gradient dT(x, 0)/^y at the wall. Under the stationary conditions 
stipulated, the heat transfer rate q w = - kw dT(x, 0)/dy is balanced by 
radiation from the surface and/or by the heat absorbed by the coolant 
as it is transported across the wall. 

As in Reference 11, the following nondiinensional similarity variables 
are introduced: 



f(n) 


tfcy) 

N / U 


V x 

oo oo 


(9) 
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The continuity equation (1) is identically satisfied by the stream 
function \|r for the two-dimensional case at hand. Equations (8) and (9) 
then determine the velocity components as follows: 


= — — fl fro 

p by ~ 2 c(t)) 


(10) 


V 



P oo 1_ / U eo v oo ' 
p 2 J x 


C f (l]) - Tlf'(Tl)], 


(ID 


where C(t]) = \Jr i|r . Equation (11) implies 


£<o) • £„ 



(11a) 


The similarity transformation given by equations (8) and (9) implies 
v w (x) = x -1 / 2 . The momentum equation (2) becomes 

f(Tj) ^ [C-^TJ) f'(Tl)] + [C-^TJ) f'(Tj)] = 0. (12) 

Additional similarity conditions are introduced by 

w = w(r]) and T = T(t]) (13) 


so that w = const, and T-. = const. Equations (3) and (4) then take 
w w 

the form 


f(Tl) 


fofri ? + jL (. . 1— 

dq dt) \S c (q) 


dw (nA 

dr] J 


0 and 


(14) 


<P C W ti 


+ i d 

r% < 1 ) d /itnA] 

k - 1 

+ 00 

"d 

V Too J Proo dT l 

CP^(T)) dT) V )_ 

+ 4 ^ 




Cpa f (T]) dw(n) _d_ / t(tiA 

S c (t]) dT l dT l \ T oo / 


= 0, 


(15) 
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where the relation a m = c Poo (K^ - 1) has been used. The coefficients 
cp^, cp , cp c , and cp k in the system of differential equations (12), (14), 
and (15) depend on w(q) and T(t}). The coefficient cp a ^ is a function of 
T(t]) only. The boundary conditions (5) and (6) take the following form: 


f (0) = 0, f (0) = f w , w(0) = w w , T(0) = T„, 


lim f' (ri) = 2, lim w(q) = w^, lim T(q) = T^, and 

T] — > 00 T] — > 00 J) oo 

f (o) y.'122 

t(U) S c (0) [1 - w(0) ] * 


(16) 

(17) 

(18) 


Equations (12) and (14) through (17) define a split boundary-value 
problem for a system of three nonlinear ordinary differential equations 
of the seventh order. Solutions with regard to Pr = Sc = 1 and 
Hx> = cp a f “ 0 are particularly simple, but constitute rather poor 
approximations to solutions with variable properties. Under the 
assumptions of = 0 and constant properties, the authors have solved 
equations (12) and (14) through (18) in a wide range of Pr, Sc, c, and 
f w . A correlation formula for the results is presented in Appendix B. 

A numerical solution of this problem for any given set of coefficients 
cp p (T, w), cp p (T, w), cp k (T, w), cp c (T, w), and cp a f ( T » w) requires an 
iterative technique which employs the four given initial conditions 
(16) and three iteratively corrected initial assumptions of f" (0) , 
w f (0), and T ! (0) to satisfy the conditions (17). The boundary values 
at the wall have to be consistent with equation (18). It is desirable 
to uncouple these differential equations in order to reduce the amount 
of numerical work involved in the solution procedure. 


III. SOLUTIONS OF THE DIFFERENTIAL EQUATIONS 

The differential equations (12), (14), and (15) can be uncoupled 
if the following relations hold true: 


c(n> i, 

*00 ^oo 

(19) 

Cp k <T)) k (n) M-co 

cPpOi) k o 0 ^(n) =1 ’ 

(20) 



CpC< T l) = 


( 21 ) 


7 


Sc^l = ^ (n) = j d 

Scoo p( n ) D(q) p* “ A ’ and 


( 22 ) 





Cp (tj) - Cp f ( Tj) 
C - C 




5 = 1 . 


(23) 


Conditions (20) and (21) imply Pr(q) = Pr^. 

Figure 10 presents the relation between f" (0) and f(0) which follows 
from numerical solutions of equation (12) under the assumption of C (rj) = 1. 
If f w is prescribed, the three initial conditions are known and thus equa- 
tion (12) can be integrated as an initial value problem. Equations (12) 
and (19) yield the relation 


f f <TJ> dTj = - In . 


(24) 


Equations (14), (22), and (24) give for Sc = const. 


w(tj) - w(0) = w' (0) 



and 


w' (0) 



(25) 


(26) 


The parameter w' (0) can be eliminated between equations (18) and (26) 
so that 


where the right-hand side is a function of f(0) = f only. Figure 11 
presents the relation between w T7 and f T7 which follows for w^ = 0 from 
equation (27) . 

Equation (15) can be written to give 


t" (tj) + t' (tj) 


Pr 

Pr f<„) - ^ 


k - 1 


Cp w' (t]) 
a f 


+ Pr 


Mf 


‘d 2 f 


dr) 


& 


2 


( 28 ; 


where t(q) = T(q)/T 00 and Pr, Sc, and cp a f are constants. The short 
notations 


Pr 

■ sJ fa, 


w,, - w, 


00 “ w w 


koo - 1 


and 


Sc 


a 2 = Pr — [f"(0)] Pr 


tlinl 

f"(0) 


dq 


(29; 


are employed, which are defined since f" (0) ^ 0. Integration of equa- 
tion (28) gives 


Pr 






t' (T|) - 


f" rri'i 


exp 
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In view of the boundary condition lim t(T}) = 1, a quadrature of equa- 
tion (30) yields Tj — >oo 



O O 


Equation (31) shows that the contribution of dissipation to heat 
transfer is governed by the parameter (k - 1) T / (T - T(0)). The 
function f" (r))/f" (0) in the right-hand side of equation (31) is obtained 
by integrating equation (12). Equation (31) includes the case of zero 
mass transfer (f w = w w = w B = 0) , m 0. Because of the boundary 
condition (19), equation (12) gives in this case for constant properties 


f"(0) = 1.328. 


(32) 


For f w = = 0, both and a 2 equal zero. Equation (31) then gives 

the well known expression 


r_i_ 

rT(ri) " 

l / Vco _ 

Nu 1.328 Pr 


_T W - Too 

r -7 L 7u ot l 

i1=n 

'•/5T 2 


1 

00 

f [f"(ri)] Pr dq 

o 


(33) 


See Chapter XIV of Reference 24. 
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In two special cases, the repeated quadratures in equation (31) 
can be reduced to one single quadrature in the expression for C^. If 
Pr = Sc = 1, 



- Q!i 4- 2 


T m 



exp 


2q-i 

f"(0) 


1 


f"(Q) ag T °° 

«l Too - T w ' 


(34) 


If Pr = Sc and M^ = 0, 


d T(rt) Top 
Tw - Too 



-l Tj=0 



- a-j 


exp 


cPa f (w °o - w w ) 



(35) 


The heat transfer rate for any values of Pr, Sc, and can be expressed 
as follows by use of equation (8): 


q w = k(o) 


dT (0) dr\ (0) = J_ k(0) f xUqq 
dr) dy 2x \±(0) ^°° V 


dT(0) 
dq • 


(36) 


The numerical values of p(0) and lc(0) may be taken from Figures 2 and 3, 
respectively, if the injection of H 2 0 vapor, He, or H 2 into air is con- 
sidered. 


IV. DISCUSSION OF THE VALIDITY OF THE ASSUMPTIONS 
(19) - (23) ON THE MIXTURE PROPERTIES 

Mixture Properties 

The simultaneous system of differential equations (12), (14), and 

(15) can be uncoupled if the binary mixture flow represented by these 
equations satisfies conditions (19) through (23) within a reasonable 
margin of error for a triple of boundary values T^, T w , and w w . Clearly, 
these conditions can be satisfied rigorously by assuming constant 
mixture properties. Condition (23) can always be met in an approximate 
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way because the specific heat c of any gas increases together with 
temperature, e.g., Figure 4. Ir > T w and w^ < w w , the conditions (19) 

through (22) can be satisfied only if each one of the mixture parameters 
p.p, p/k, Cp, and Sc is influenced in the same direction by changes of 
either temperature T or concentration w. Results of pertinent methods 
listed in Appendix A show that this situation prevails if H 2 0 vapor, He, 
or H 2 is injected into air (see Figures 5 through 9)* This condition of 
the same influences exerted by T and w fails to be satisfied for several 
investigated foreign gases with molecular weights exceeding the one of 
air. 


As a first approximation, conditions (19) through (23) are relaxed 
and satisfied by equating values at the wall and at the outer edge of 
the boundary layer. By employing the T- and w-dependencies of the 
mixture properties presented in Figures 1 through 5, these relaxed 
conditions have been satisifed approximately for given temperatures 
and T w listed in Table I by calculating suitable values of w w . Equa- 
tions (25), (26), (30), and (31) have been evaluated for T(q) and w(tj) 
with regard to case 3 listed in Table I (see Figures 12 and 13). The 
functions T(tj) for Moo = 0 and w(q) presented in Figures 12 and 13, 
respectively, have been used to evaluate conditions (19) through (23) 
in the range 0 - q < oo. The results of this calculation appear in 
Figure 14. At least in this one special case, the margin of error in 
satisfying conditions (19) through (22) is smaller inside the boundary 
layer than at the wall. It may be conjectured for M^, = 0 that it is 
sufficient to satisfy the relaxed conditions (19) through (23) by employ- 
ing only given boundary values at the wall and at the outer edge of the 
boundary layer, instead of the unknown functions T(t]) and w(rj) . A glance 
at Figure 13 shows that the relaxed form of conditions (19) through (23) 
is not applicable unless the difference between the temperature profiles 
for a given Mach number M^ and for M^ = 0 is sufficiently small. This 
consideration is verified by the final results presented in Figure 15. 


V. RELATIONS FROM THE REFERENCE ENTHALPY METHOD USED FOR COMPARISON 

The reference enthalpy method derived by Eckert in Reference 7 
gives the expression 


%jo ^i 


<*r - 


i w> 


(37) 


for the heat transfer at an impermeable solid wall for compressible air 
boundary layer flow with constant pressure. The recovery enthalpy can 
be expressed as follows 


i 

r 



O^oo 


1) Too. 


(38) 
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For laminar air flow, the recovery factor can be approximated by 


r . = s/Pr . 


( 39 ) 


The temperature-dependent gas properties are introduced at a reference 
enthalpy defines as 


i* = i + 0.5 (i - i ) + 0.22 (i - i ) ; (40) 

i* gives the reference temperature T* by use of enthalpy tables for air. 
For two-dimensional laminar flow along a surface with constant properties, 
pressure, and temperature, the local heat transfer coefficient is obtained 
from the relation 


h. 

i 


P u „ 


0.332 

Re l/2 p r 2/ 3 ' 


(41) 


It has been shown in Reference 7 that equation (41) is a correlation 
formula for heat transfer in high-speed flow with variable properties if 
p, Re, and Pr are evaluated at the reference temperature T*, i.e., 

Re (x) = p'^xu^/ij* and Pr = pi*c p */k*. Equations (36) through (41) give 
for Wo, = 0 


k(T w , wj /p(Tj n(T 00 ) r p(T x ) c n 0T) 


2/3 


w 




(T , w„) 


w 


w y ± p(T“ ) u(T“) 


g. 

k(T-) 


dT(Q) 

dq 


%o 


0.664 


|j(r) c p (rQ M; 


2 C P, 


(kco - 1) T ro - i^ 


(42) 


For constant pressure, it has been shown in Reference 14 that 
exact results for heat transfer in the presence of mass transfer cooling 
with a foreign gas (e.g., References 2, 8, and 26) can be approximated 
by the following expression: 


^wo 


1 


- 1.82 



p v r p xu t ~ 

’w W / 00 oo 1 

Pco U oo J M* C* 


(43) 


13 


is given by 


c * . j a fflQ. h <t*) 

p. 

Equation (43) is valid over a wide range of values M^, T^, T w , and P w v w . 
Figures 20 and 21 in Reference 14 indicate for the injection of H 2 and 
He, respectively, that evaluations of equation (43) deviate by less 
than approximately 0.05 from rigorous results derived by other investi- 
gators. These rigorous results and pertinent experimental data show, 
however, that q™/ q wo ceases to he a linear function of the mass transfer 
rate p w v w as q w /q WQ 0. For this reason, it is questionable whether or 
not equation (43) is valid for large mass transfer rates pwVw* 


VI. RESULTS OF THE DERIVED METHOD 

Table 1 presents triples of numbers Too, w w , and T w for which the 
relaxed form of conditions (19) through (23) has been satisfied within 
the margin of error also listed in this table. The slopes of both pja 
and p/k as functions of w in Figures 6 and 7, respectively, show why 
the cases appearing in Table 1 cover a wide w w -range for H 2 0-air mixtures 
and only small values of for He-air and H 2 -air mixtures. The width 
of the bands presenting Sc and Cp as functions of w for T - const, in 
Figures 9 and 5, respectively, confirms this situation. 

When the five conditions (19) through (23) have been satisfied in 
the outlined approximate way, suitable average values of Pr, Sc, and 
cp a £ can be estimated. For each one of the cases listed in Table I, 

Figure 11 has been used to determine f w as a function of Sc and w w . 

Figure 10 gives f^ in terms of f w . The momentum equation (12) can then 
be integrated as an initial value problem to give f(q). The resulting 
function f" (q)/f" (0) is substituted into equation (31) to obtain T 1 (0) 
by quadratures. Equation (42) finally yields q w /q wo - For = the 
numerical evaluations involved in this method can be reduced to a few 
hand calculations by employing the correlation formula (B-2) in 
Appendix B. 

The marked points in Figure 15 present results of equation (42), 
which is a consequence of the derived method, for the H 2 0-air mixtures 
listed in Table I. The curves in Figure 15 have been obtained from 
the correlation method presented in Reference 14. Equation (43) gives 
the upper line when evaluated for the molecular weight ratio of H 2 0 and 
air. The lower curve represents a special correlation formula shown in 
Figure 22 of Reference 14 for H 2 0 and air. The agreement between points 
referring to M^ = 0 and the lines is remarkably close, particularly in 
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view of the accuracy of the correlation method of Reference 14, where 
rigorous solutions for H 2 0 and air cover a band of certain lateral 
width around the lower straight line, which is also approximated by 
the upper straight line representing equation (43) within a certain 
margin of error. While the agreement between the derived method and 
the correlation method of Reference 14 is good for H 2 0-air mixtures 
in a low-speed boundary layer, the evaluations of equation (42) differ 
significantly from the lower straight line for = 3. This is due to 
the fact that the influences of T(r|) and w(q) on the pertinent property 
parameters given by equations (19) through (23) do not compensate each 
other in an approximate way as they do for = 0, since the temperature 
profiles reach a maximum inside the boundary layer for = 3 (see 
Figure 13). This trend is particularly strong when Tw and Too differ 
only by a relatively small quantity, see case no. 2 in Figures 13 and 
15. Points for considerably higher mass transfer rates f w could have 
been calculated by specifying larger values of (T^ - T w ) , This is due 
to the fact that the values of w w which satisfy conditions (19) through 
(23) increase together with (T OT - T w ) . 

VII. REVIEW AND APPLICATIONS OF THE DERIVED METHOD 

A laminar air boundary layer with injection of foreign gas is con- 
sidered here. The mixture properties may depend in any prescribed way 
on temperature T and foreign gas concentration w. The partial differential 
equations representing conservation of matter, injected species, momentum, 
and energy can be reduced to three simultaneous ordinary differential 
equations by a similarity transformation which^ is valid both at the stag- 
nation point and for a constant pressure wall. Several investigators have 
solved this simultaneous system of similarity equations by rigorous numeri- 
cal methods for the case of constant pressure. The solution procedures 
employed are quite cumbersome because a triple iteration is necessary to 
satisfy three boundary conditions at the outer edge of the boundary layer. 
Available rigorous numerical solutions have been correlated in Reference 14 
and represented by a linear dependency between q w / q wo and the injection 
rate / p^u^, whose slope coefficient is a function of M a /M£, Re (x) , 

and pV/Pbfe. While this correlation formula treats the case of con- 
stant pressure in a satisfactory manner, only individual numerical solu- 
tions for special given values of T^, M^, and T w have been published so 
far for the stagnation point. 

Because of this situation, a relatively simple calculation method 
is proposed in this paper which gives q w / q wo for both the stagnation point 
and the case of constant pressure. To show its performance and limita- 
tions, this method is worked out here only for the latter case. The 
method rests on the fact that the three ordinary differential equations 
in similarity variables can be uncoupled if the mixture parameters pip, 

pi/ k. Sc, and c p , and (c - c ) are constant. To render these parameters 
F F a Ff 
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independent of the solutions T(q) and w(q) to be determined, these five 
conditions are relaxed so that only the values at the wall are equated 
to the values at the outer edge of the boundary layer. Comparison to 
the correlation solution compiled in Reference 14 shows that an approxi- 
mate matching of the relaxed form of the five conditions is sufficient 
unless dissipation has a significant influence on the temperature pro- 
file, These relaxed conditions can be satisfied in an approximate way 
for the injection of H 2 0, He, or H 2 into air because changes of either 
T or w have the same influences on each one of the mixture parameters 
[ip, p/k, Sc, and c p in these three cases. This situation does not 
exist, however, for several investigated foreign gases whose molecular 
weight exceeds the one of air. The relaxed form of the conditions can 
be satisfied in the approximate manner for He or H 2 if the injection 
rates are very small. In case of H 2 0 injection into air, however, 
values of q w / q wo as low as 0.5 can be obtained for moderate values of 
* T w ). If this temperature difference is sufficiently large, the 
mass transfer rate p w v w / p^u^ of H 2 0 vapor may reach values outside the 
range of applicability of the correlation formula (43) which is linear 
in pw^w/ Pco^oo • 

If the five conditions on the mixture properties are satisfied in 
the outlined approximate way, the heat transfer parameter q w / q WQ is 
obtained by integration of equation (12) and quadratures in the right- 
hand side of equation (31). For - 0, the numerical evaluations can 
be reduced to a few hand calculations by employing the correlation 
functions presented in Appendix B for the Euler numbers m - 0 (con- 
stant pressure), m = 0.5, and m = 1 (stagnation point). 

In conclusion it may be pointed out that the method derived in 
this paper can give simple engineering solutions of the heat and mass 
transfer problem which in some applications exceed the range of validity 
of the correlation formula (43) of Reference 14. Evidently, the derived 
method is valid for primary gases other than air and for foreign gas 
concentrations w^ # 0. The method can be generalized to include homo- 
geneous or heterogeneous chemical reactions of given order. Also, the 
method can be rewritten to account for a mass transfer vector component 
tangential to the wall. 
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APPENDIX A 

THE MIXTURE PROPERTIES 

1 . Viscosity, Diffusion, and Thermal Conductivity Coefficient 

Viscosity, ordinary diffusion, and thermal conductivity coef- 
ficients were determined by use of a rigorous kinetic theory of gases 
of low density initially developed for pure gases by Chapman (Refer- 
ence 6) and independently by Enskog at the same time. Hirschfelder and 
associates (Reference 19, pp. 514-540) extended this theory to include ^ 
binary mixtures. The Chapman-Enskog theory expresses the transport 
coefficients in terms of the potential energy of attraction between ^ l 
a pair of gas molecules. Since exact potential functions ar e/lcnofalf^ ' 
empirical relations are commonly used to describe interaction^ between 
molecules. According to Reference 19, transport properties of pure 
gases can be expressed in terms of molecular weight, pressure, tempera- 
ture, specific heat, a characteristic diameter a of the molecule, and 
collision integrals which represent the deviation of the molecules from 
the rigid sphere model. The collision integral is a function of a 
dimensionless temperature (kT/g) and depends upon the potential model 
under consideration. The parameters e/ K (°K) and o(A) are often referred 
to as "force constants"; c is a characteristic energy of attraction 
between two molecules, and K is the Boltzmann constant. 

In calculations of the viscosity and thermal conductivity of 
H 2 and He, tabulated values of the collision integrals based on the 
Lennard- Jones potential model were taken from Reference 19, pp. 1126- 
1127. Force constants were also obtained from Reference 19, p. 110. 

The expression used for calculating thermal conductivity was developed 
on a semi-empirical basis by Eucken (Reference 12) and was also derived 
theoretically by Hirschfelder, Reference 19, p. 534. 

Hirschfelder 1 s general mixture equations of viscosity, 

Reference 19, p. 531, mass diffusivity, Reference 19, p. 539, and 
thermal conductivity, Reference 19, p. 535, were used in computations 
for H 2 -air and He-air mixtures. No additional input data are needed 
for the use of these equations except the mole fraction of foreign gas 
and the force constants of the mixture. Empirical combining laws were 
used to obtain mixture force constants which account for interaction 
between unlike molecules, 


°12 = 0-5(0! + a 2 ) and e 12 = nTcj 


(A . 1 ) 


(see p. 222 of Reference 19). 
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Approximate, semi- empirical relations were used as a comparison 
method for determining viscosity and thermal conductivity of H 2 -air and 
He-air mixtures. The expression used for viscosity calculations (Refer- 
ence 5, p. 5) was first developed by Buddenberg and Wilke (Reference 4, 
p. 1345) and simplified according to kinetic theory by Wilke (Refer- 
ence 28, p. 57). In order to evaluate Wilke's equations, the tempera- 
ture, concentration, molecular weights, and viscosities of the pure 
gases are needed. The pure component viscosities were obtained by use 
of Hirschfelder' s modification of the Chapman-Enskog theory outlined 
previously. A semi-empirical method for thermal conductivity derived 
by Enskog and presented by Carlson and Schneider (Reference 5, p. 13) 
was used in the same manner as the Wilke method for viscosity. Perti- 
nent tables in Reference 10 have been used for comparison. 

Computations for the viscosity and thermal conductivity of water 
vapor were based on the Stockmayer potential function which is a more 
realistic model for polar molecules than is the Lennard- Jones model. 

The collision integrals for this case were taken from tables in Refer- 
ence 22, pp. 23-33. Force constants for water molecules were obtained 
from Reference 19, p. 599. Force constants for H 2 0-air mixtures were 
determined from equation (A-l) and were then multiplied by a correction 
factor (Reference 19, p. 600) to account for the effective potential 
energy of attraction between a polar and nonpolar molecule. The 
viscosity and thermal conductivity for H 2 0-air mixtures were then 
determined by the theoretical method used for H 2 -air and He-air mix- 
tures. Diffusion coefficients for H 2 0-air mixtures were calculated 
according to an empirical equation given by Eckert (Reference 9, p. 512). 
Figures 2, 3, and 8 present viscosity, thermal conductivity, and 
ordinary diffusion coefficient, respectively, as functions of foreign 
gas concentration w for selected temperatures T. 

2. Specific Heat and Density 

Experimental values of the specific heats of pure gases 
(References 13, 18, and 20) were used in all cases except in calcula- 
tions of the thermal conductivity of H 2 -air and He-air mixtures (see 
Figure 4). For this case, the specific heats of the pure components 
were determined by use of the energy principle of equipartition. The 
actual contributions of the various energies for each degree of freedom 
of the molecule to the total energy of the molecule were estimated accord- 
ing to expressions given by Hirschfelder in Reference 19, pp. 116-119. 

The specific heats of mixtures were determined by the following ideal 
gas relationship: 


c p (T, w) = Cp^(T) w + C p (T) (1 - w). 


(A- 2) 


I 
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(Reference 5, p. 3). Densities of pure gases and mixtures were deter- 
mined from the following ideal gas relationship: 


Pi = 


p.M. 

li 

RT 


and 


PM- 


RT[1 




«< 5 ? - 


(A- 3) 


(Reference 5, p. 3). The average molecular weight of air was taken 
as 28.97 (79% N 2 and 21% 0 2 ) . Figures 1 and 5 present the density 
and the specific heat as functions of foreign gas concentration w for 
selected values of temperature T. 
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APPENDIX B 

CORRELATION FORMULAS FOR NuA/Re IN CASE OF Mm = 0 

The similarity equations (12) and (14) through (18) have been 
solved numerically on a digital computer under conditions of Mm = 0, 
constant properties, and Euler numbers m = 0, m = 0.5, and m = 1. The 
Euler number is defined by u^x) = c x®; m = 0 and 1, respectively, 
represent the cases under discussion in this paper, i.e., constant 
pressure and stagnation point. The numerical results for the heat 
transfer parameter 


Nu/s/Ri = - x °^ /(T w - TjJte 


have been correlated as functions of the pertinent constant parameters, 
Pr, Sc, c = Cpf/Cpa, and B where B = - f w Pr(m + l)/2 [NuA/Re]^ is 


the blowing factor. The correlation functions are valid in the follow- 
ing bounds: 



g 

g 

g 

£ 


B 

Pr 

c 

Sc 


g 

g 

g 

< 


0.1 

5.0 

5.0 

5.0 


(B-l) 


For 0.25 § c § 1.0, the blowing factor may vary in the bounds 
0 g B S 0.3. For m = 0, 


NuA/Re = 0.33203Pr°* 3489 - Bg^Sc, c)Fi(Pr) -jl + S^Pr, c) H x (c) + 


+ A fScl A (c) - £ 1 (Pr , c ) 

i (S ^ Az(C) g(Sc, c) F x (Pr) 


(B-2) 


where 


F x (Pr) = 0.835 + 0.18 Pr, 
g 1 (Sc, c) = 0.2394 + (c - 1) 


(-0.00312 c + 0.18213 + Sc (0.00059 c 


- 0.03298) + Sc 2 (-0. 00005 c + 0.00364) 
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Si(Pr, 

»i(c) > 
Ai(Sc) 
A 2 (c) = 
Ei(Pr, 
For m = 
Nu /\l Re 

where 

F 2 (Pr) 

g 2 (Sc, 

H 2 (c) = 

G(c) = 

E 2 (Pr) 

For m = 

Nu /si Re 


N (Pr - 1) (c - 1) 

c ) - - | Pr _ i I I c - 1 I and 0 if Pr = 1 or c = 1, 


= c - 1 


0.115 
0.63 + c 


= Sc 


0.21066 - 0.05479 c + 0.01058 c 2 , and 


. 0.6407 

c) = c -0.33 + 2.19 c -!- 0.63 Pr + 0.25/Pr 


10 


-3 


= 0.5, 

0 3777 f 

= 0.4775Pr * - Bg 2 (Sc, c) F 2 (Pr) jl + H 2 (c) + G(c) E 2 (Pr) 


(B-3) 


= 0.76 + 0.265 Pr, 
c) = 0.3036 + (c - 1) 


0.25633 - 0.05076 Sc + 0.00583 Sc 2 ] 

J 


: (1 - c) 


0.27902 - 0.44556 c + 0.18848 c 2 - 0.02194 c 3 


(c - 1 ) / | c - 1 and 0 if c = 1, and 


* 0.40581 - 0.29358/Pr - 0.05363 Pr + 0.02649/Pr 2 . 


: 1 . 0 , 


0.57047Pr°* 3858 - Bg 3 (Sc, c) F 3 (Pr) ^1 + S x (Pr, c) 


+ E 3 (Pr) 


H 3 (c) + 


(B-4) 
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where 

F 3 (Pr) = 0.762 + 0.275 Pr, 


g 3 (Sc, c) = 0.3577 + (c - 1) 


0.30035 - 0.05438 Sc + 0.00589 Sc 2 , 


H 3 (c) = (1 


c) 


0.26302 - 0.30873 c + 0.12054 c 2 


0.01369 c 3 , 


and 


E 3 (Pr) = 0.45561 - 0.32596/Pr - 0.06026 Pr + 0.02916 Pr 2 . 


Comparison of the computer solutions of equations (12) and (14) through 
(18) to the correlation formulas (B-2) through (B-4) shows a maximum 
error of 0.032 of NuA/lle in the bounds def ined by the relations (B-l). 
NuA/Re varies in the bounds 0.158 ^ Nu/s/Re ^ 1.025 in this region. 

The average error of Nu/s/Re in this region is approximately equal to 

0 . 01 . 


NOTE: The correlation formulas presented in Appendix B have been 

derived by Mr. F. E. McKinney of the George C. Marshall Space 
Flight Center. 
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FIG. 12. CONCENTRATION PROFILES 
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FIG. 13. TEMPERATURE PROFILES 
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